Air expansion in the water rocket 
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We study the thermodynamics of the water rocket in the thrust phase, taking into account the 
expansion of the air with water vapor, vapor condensation and the energy taken from the environ- 
ment. We set up a simple experimental device with a stationary bottle and verified that the gas 
expansion in the bottle is well approximated by a polytropic process PV I3 = constant, where the 
parameter /3 depends on the initial conditions. We find an analytical expression for f3 that only 
depends on the thermodynamic initial conditions and is in good agreement with the experimental 
results. 
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I. INTRODUCTION 

The water rocket is a popular toy that is already sev- 
eral decades old Its inexpensive version can be built 
with a plastic bottle filled with water and pressurized 
air obtained with a bicycle pump. When the bottle is 
opened, the internal air pressure pushes the water out 
and this causes an increase in the momentum of the bot- 
tle. The water ejection gives the bottle the thrust that 
allows it to move for several meters, see Ref. [2j. 

The launch of water rockets is used in our undergradu- 
ate laboratory to illustrate several physics concepts, such 
as : (a) Newton laws, (b) conservation of momentum, (c) 
work-energy theorem, (d) Bernoulli equation, (e) mass 
conservation equation and (f) ideal gas expansion. On 
this last point it is necessary to clarify whether the ex- 
pansion of the gas in the rocket is an adiabatic or an 
isothermal process. 

The theoretical description of the water rocket is not 
trivial and in general it is necessary to adopt certain 
premises to obtain a friendly caricature of reality. In 
particular, the rocket's only source of energy, the air ex- 
pansion, is in general modeled 043 as an isothermal or 
an adiabatic process involving dry air. The qualitative 
predictions of these theoretical approaches give reason- 
able results, however their quantitative description are 
incorrect, showing that some hypothesis of the model is 
inadequate. In particular observing the launch of the 
water rocket, we see that when the ejection of the water 
concludes, an additional ejection of fog follows. There- 
fore, it is clear that the expansion occurs for a mixture of 
dry aii\ water vapor and condensed water (fog). In refer- 
ence p the author deepens in this point and he builds a 
pressure- volume relation assuming that the total entropy 
of the mixture is conserved during the expansion; he ad- 
ditionally shows that the solution of the pressure- volume 
equation can be approximated by a polytropic process. 



We developed an experimental device in order to elu- 
cidate the nature of the water-rocket gas expansion. We 
also present, a brief theoretical development that pro- 
vides an analytical expression for the parameter that 
characterize the polytropic processes. 

In the following we shall use the term "gas" to refer 
to the mixture of: dry air and water vapor. The paper 
is organized as follows, in the next section we study the 
exchange of heat between the gas and the environment 
in a polytropic process. The experimental device and 
our results are presented in the third section. In the last 
section some conclusions are drafted. 



II. POLYTROPIC PROCESS AND EXCHANGE 
OF ENERGY 

We present here a brief revision of the polytropic pro- 
cess and the calculation of some quantities of heat in- 
volved. A polytropic process is a quasistatic process car- 
ried out in such a way that the specific heat remains 
constant during the entire process [7|-[9|. Therefore the 
heat exchange when the temperature changes is 



dQ 
dT 



Nc, 



(1) 
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where Q is the heat absorbed by the system, T is the tem- 
perature, c is a constant specific heat and N is the num- 
ber of moles. Pressure and volume can change arbitrarily 
in this process, therefore the specific heat also depends of 
the process. Of course, the processes in which the pres- 
sure or the volume are kept constant are polytropic with 
specific heats c p or c v respectively. In an adiabatic evo- 
lution there is no heat exchange between the system and 
the environment, therefore it is a polytropic process with 
c = 0. At the other extreme, the isothermal evolution is 
characterized by a constant temperature therefore it may 
be thought as a polytropic process with infinite specific 
heat. 

Using the first law of thermodynamics together with 
the internal energy of an ideal gas U = c v NT and Eq.(TT]), 
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we have for an infinitesimal polytropic process 

(c v - c)NdT + PdV = 0, (2) 

where P is the pressure and V the volume. From Eq.([5]) 
and using the ideal gas equation of state, it is easy to 
prove that the polytropic process satisfies 



where 



PV^ = constant, 



(3 = ^ 



(3) 



(4) 



Usually thermodynamic textbooks define a polytropic 
process through Eq.(j3]) (for example [13]). However we 
have started with Eq.(UJ which gives a direct link be- 
tween heat and temperature changes. In particular for 
an adiabatic process c = and j3 — 7 = c p /c v and for 
an isothermal process, c = 00 and j3 = 1. For the pro- 
cess treated in this paper we find that 1 < j3 < 7 where 
7 = 1.4 for dry air at room temperature. The specific 
heat as a function of /3 is 



f3c v 



1-/3 



(5) 



The total heat Q exchanged between the gas and the 
environment can be calculated integrating Eq.([T]) 



Q = cN(T - To) 



(6) 



where T is the gas temperature and To is the environ- 
ment temperature, which is also the water temperature 
and the initial temperature of the gas in the bottle. An 
experimental measurement of the time dependence of T 
determines the time dependence of Q through Eq.©. 

As the water is expelled from the bottle the gas ex- 
pands and cools. In this process the relative humidity 
increases up to saturation and some of the vapor con- 
denses. The quantity of condensed water vapor at the 
end of the experiment will depend on the initial relative 
humidity in the gas and the initial and final tempera- 
tures. We want to estimate the mass and the heat in- 
volved in this condensation. The mass of water vapor 
m v can be calculated using [l2| 



PyV 

RyT' 



(7) 



where P v is the partial pressure of the vapor at T and 
R v = 0.4615 kJ/(kg K) is the ideal gas constant for the 
vapor. The vapor partial pressure is determined by de 
equation P v — ee w where e £ [0, 1] is the fraction of rel- 
ative humidity and e w is the saturation pressure which 
depends only on the temperature and may be obtained 
from a vapor table. The pressurized air is generated by 
the compressor at a higher temperature than the envi- 
ronment and its cooling results in the saturation of the 
air with vapor inside the system. Thus we assume that 



we are working with e w = 1 during all the course of the 
experiments, then P v = e w . This assumption implies the 
maximization of the energy delivered by the vapor-liquid 
phase transition. In our case the initial temperature is 
To = 292 K and partial pressure is P VQ — 2.3 kPa (this 
last value was taken from a vapor table). 

In order to get a fuller understanding about the heat 
exchange due to vapor condensation we start with the 
Clapeyron-Clausius equation [ll[ for the vapor-liquid 
phase transition 



dp 
~dT 



R v T 2: 



(8) 



where L v is the latent heat, which for our purposes may 
be considered constant, L v = 2500 kJ/kg. The solution 
of Eq. (jHJ gives the partial pressure P v function 
of the temperature T 



tJ--±) 



P v = P Vo e^' T o 
where the constant was adjusted so that P v (Tq) = P, 



(9) 



The condensation produces a change in the mass of 
vapor in the gas. The mass of condensed vapor can be 
calculated using Eq.([7]) twice, that is 



M = 



Pv V P V V 



RyTo r v t 



(10) 



It is easy to prove that the polytropic process is also 

g 

characterized by the equation PT l -$ = constant. Then, 



V 
T 



NR 
~P~ 



NR /Tb 
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where R is the ideal gas constant. To obtain the latent 
heat of condensation Q v = ML V we substitute Eqs. (|9lll| ) 
into Eq.([l0]) that is 



NRP V0 L V 
RvPo 



i„ (i-lo) ( To 
T 



(12) 



One may wonder if there are other types of heat ex- 
change between the system and the environment. As 
will be seen our experiment lasts for less than one sec- 
ond, then it is very reasonable to assume that the heat 
exchanged through the walls is not important, the energy 
delivered by the vapor-liquid phase transition remains as 
the sole source of heat for the system. 

It is important to point out that, due to the water- 
vapor condensation, the number of moles N in the gas is 
not constant. However, Eqs. (l6ll2l) were obtained using 
the ideal gas law for air and water vapor, which assume 
that N and N v (the number of moles of water vapor) 
are constants. Therefore, Eqs. (|6ll2j) may be used only if 
the variation of moles is negligible during condensation. 
The relative number of moles condensed 5N V /N V can be 
estimated using Eqs. (|7ll0[) 



SN V 
N v 



M 



= 1-1^ 
T 



4» ( 1 . 
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(13) 
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where m Vo is the initial water vapor mass. If we esti- 
mated the final temperature as T = T — 10 K then 
Eq.(dHJ gives us the value 6N V /N V ~ 0.6. The total 
relative variation of moles in the gas SN V /N is 

<5iV„ 8N V N v 



Pv SN V 



N 



N v N P N v 



10- 



(14) 



where we have used the ideal equation gas both for the 
gas as for the water vapor. Therefore we have the fol- 
lowing situation, SN V /N is negligible but SN V /N V is not. 
In the deduction of Eo. (fT2)) . we have used Eq.© which 
is correct only if SN V /N V ~ 0. Then, it is clear that 
Eq.(fT2j) can only be correct in the limit T — > Tq. There- 
fore, we must substitute Eq. dT^t by the following exact 
differential equation 



, n _ L v P V0 NR 



P 



T R V P \l—/3 R v T 



dT, 



(15) 



that it is obtained replacing T = Tq + dT in Eq. (JT2J) and 
taking the terms to first order in dT. Now it is possible 
to equate dQ v with dQ taken from Eq.([T]) 



dQ = dQ v , 
and using Eqs. (|5|15[ ) in Eq. (fT6|) lead to 

Cp/R- f3c v /R _ L v P Vg f f3 L v 



1-/3 



T R V P \ 1 — p R v Tq 



(16) 



(17) 



This last equation gives a theoretical expression to obtain 
the constant parameter /3, 



P 



c p /R+(L v /R v T y(P Vo /P ) 



c v /R + (L v /R v T y (P V0 /P ) - (L V /R V T ) (P V JP ) 

(18) 

This is the main theoretical result in this paper and from 
its analysis we conclude: (a) If P Vo ~ or L v ~ then 
P ~ 7 = this means that the process is adiabatic 
for dry air. (b) ft depends only on the thermodynamical 
initial conditions of the gas, this is unexpected since other 
parameters like the water ejection time do not appear. 
Tab.Q]shows the theoretical values of (3 Eq. ([T5| . which we 
call f3* , computed for our experimental initial conditions. 

It is also interesting to point out that the model devel- 
oped in Ref. Q also predicts that the exponent (3 depends 
mostly on the initial relative humidity and proposes the 
empirical function 



jfl = 1.15 + (1.4 - 1.15) cxp(-36P„ /fb). 



(19) 



Fig.flT]) shows P (Eas. (|18ll9p ) as a function of the initial 
ratio P Vo /Po (that is the molar fraction), the figure also 
shows the experimental values obtained in the present 
paper. 




FIG. 1: The polytropic exponent as a function of the initial 
molar fraction. The star points are the experimental values 
of 13. The solid line is our theoretical calculation given by 
Eq. (|18|) and in dashed line is the function given by Eq. (|19l) . 



device to measure the heat exchange between the com- 
pressed gas of the water rocket and the environment. The 
water rocket is fixed to the laboratory frame and then the 
time dependence of the pressure and volume of the ex- 
panding gas in the bottle is measured. The stationary 
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FIG. 2: The stationary bottle, (A) is a tire valve, (B) is a 
III. EXPERIMENTAL RESULTS pressure manometer, (C) is a faucet, Z is the height of the 

gas column situated over the water column. 

In order to elucidate the nature of the expansion of 
the air in the water rocket we designed an experimental bottle is depicted in Fig. [21 The bottle is a transparent 
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acrylic cylinder of inner diameter 5.40 cm, wall thickness 
0.30 cm and length 1 m, that is maintained in vertical po- 
sition. The upper and lower covers are both of aluminum. 
The upper one has a tire valve and a manometer (WIKA 
EN 837 - 1) working from to 4 bar. The lower cover 
has a conical shape with a slope of 30° with a 0.5" faucet 
that allows to fill and evacuate the water in the bottle. 
The bottle is was built in our mechanical workshop. The 
device is quite simple and the experiment is not danger- 
ous for the range of pressures used in this paper, so it 
is appropriate for undergraduate students. In a typical 
trial the bottle is filled with about 1400 cm 3 of water and 
the rest with air. It is placed in vertical position with the 
faucet closed. The compressed air is introduced through 
the valve and then the faucet is opened. The level and 
pressure of the gas is filmed with a PIXELINK PL-B776F 
camera. The video was obtained at 50 frames per second, 
with a exposure time of 10 ms per frame. The ejection 
of water from the bottle lasted between 0.64 s to 0.88 s, 
depending of the initial pressure, see Tab. Q] Finally the 
filming was processed manually to obtain the pressure 
and volume of the gas as functions of time. 

Fig. [3] shows the results for five sets of experimental 
data of pressure and volume. Each set corresponds to 
different initial pressures with the same initial volume; 
for each initial pressure the experiment was repeated five 
times and all the data are incorporated in the same figure. 
The data are adjusted with a power law PV& ^constant. 
The values of /3, the slope of the adjustment in the fig- 
ure, are given in Tab. |TJ The experimental values of the 
polytropic exponents (5 are in very good agreement with 
their theoretical values j3* given by Ea. (|T8|) . 




FIG. 3: The dimensionless pressure as a function of the di- 
mensionless volume of the gas in a log-log scale. The experi- 
mental data for five experiments are shown together with their 
adjustment with thin lines. The values of the initial dimen- 
sionless pressure are, from top to bottom, 4.8,4.4,4.0,3.6,3.2, 
with Pa the atmospheric pressure and Vq = 928.6 cm 3 . 



Fa /Pa 


P 


/r 


c/R 


t(s) 


Q (J) 


4.8 


1.22 


1.27 


-2.0 


0.64 


161.7 


4.4 


1.21 


1.26 


-2.2 


0.66 


158.3 


4.0 


1.20 


1.25 


-2.6 


0.75 


154.7 


3.6 


1.19 


1.24 


-2.7 


0.75 


144.8 


3.2 


1.17 


1.23 


-3.5 


0.88 


144.9 



TABLE I: 



The fifth column of Tab|T] shows the ejection time at 
which end the water is exhausted from the bottle. The 
heat Q at the end of the process shown in the last column 
of Tab|H it is calculated using Eq.© the experimental 
data and the ideal gas equation. 



IV. DISCUSSION AND CONCLUSION 

The typical ejection time of the water in our experi- 
ments is 0.7 s and the typical ejection time in the usual 
water rocket is 0.1 s 0. The difference in the ejec- 
tion rates must be analyzed: (i) The main geometrical 
difference between the bottles is related to the nozzles. 
The nozzle section of the stationary bottle a\ and the 
standard soda bottle a<z — 3.46 10~ 4 m 2 are linked by 
a i = 0.4a2. A smaller nozzle was necessary in order to 
accommodate the capability of our video camera, (ii) 
Another difference between the experiments is the ap- 
parent gravity. In the stationary bottle the water is only 
subjected to the weight. In the water rocket, the water 
mass is subjected to a force which is the sum of its weight 
plus a fictitious force due to the rocket acceleration. The 
apparent gravitational field is several times larger than 
the Earth's gravitational field. (Hi) Finally, our theo- 
retical calculation shows that (3* does not depend on the 
water output rate, it only depends on the thermodynamic 
initial conditions of the gas. Additionally (3* agrees with 
our experimental values. Then we conclude that same 
process should occur in the water rocket with the same 
thermodynamic initial conditions. 

In summary, we have presented an experimental study 
of the air expansion in the water rocket for the first 
time, using a stationary bottle. It is shown that the 
air expansion follows a polytropic process of the type 
PV@ — constant where (3 is an average value with slight 
dependence on the initial conditions. Furthermore we 
obtained an analytical expression for the polytropic ex- 
ponent which agrees reasonably with the experimental 
values. This theoretical exponent only depends on the 
thermodynamic initial conditions of the gas and it is also 
suitable to study the air expansion in the usual water 
rocket. Finally we conclude that the vapor condensation 
is the main energy source for the gas expansion process 
in the water rocket. 
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